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The radiative evolution of the relative orientations of the spin and orbital angular momentum 
vectors Si, S 2 and L, characterizing a binary system on eccentric orbit is studied up to the second 
post-Newtonian order. As an intermediate result, all Burke-Thorne type instantaneous radiative 
changes in the spins are shown to average out over a radial period. It is proved that spin-orbit and 
spin-spin terms contribute to the radiative angular evolution equations, while Newtonian, first and 
second post- Newtonian terms together with the leading order tail terms do not. In complement 
to the spin-orbit contribution, given earlier, the spin-spin contribution is computed and split into 
two-body and self-interaction parts. The latter provide the second post-Newtonian order corrections 
to the 3/2 order Lense-Thirring description. 



I. INTRODUCTION 



Ongoing theoretical studies of the gravitational radiation emitted by compact binaries and of the radiation reaction 
on the orbit are strongly motivated by the desire to detect gravitational waves using the LIGO / VIRGO § 
interferometers. In order to extract the signal from the experimental data one needs templates of sufficient accuracy. 
Based on the expectation that during the adiabatic regime the radiation reaction will circularize the orbit || , templates 
for circular orbits are computed. If the orbits are eccentric, the detectability of the gravitational waves by using circular 
templates decreased Q . Quinlan and Shapiro [5j and Hills and Bender || argue for a significant number of eccentric 
£f) I binaries in galactic nuclei for which the time before plunging is insufficient for circularization. A faithful description of 
. the behavior of these binaries in the adiabatic regime requires not only high post-Newtonian orders, but also eccentric 
orbits. Such a generic treatment, valid up to the second post-Newtonian order was provided by Gopakumar and Iyer 

If the spins of the neutron stars /black holes which form the binary are equally considered, the number of kinematical 
variables increases significantly and the dynamics complicates. Currently a second post-Newtonian order accurate 
description of the motion of the spinning binary 0- pi| is available. Spin-orbit and spin-spin type contributions 
appear at the 3/2 and second post-Newtonian orders, respectively. There is one notable exception over this rule: the 
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generate first and 3/2 post-Newtonian o rder changes in the angles between the angular momenta vectors []. Conse- 



quences of the spin precession equations (1.1) and of the conservation of the total angular momentum J = L + Sx +S2 



holding in the second post-Newtonian order are summarized in Table 1. In the present paper the angular momenta, 



1 We evaluate the post-Newtonian order of a correction Sf to some quantity / through the following prescription. Count the 
powers of e ~ Gm/rc 2 w v 2 /c 2 , representing one post-Newtonian order, in both Sf and /. (G is the gravitational constant, 
c the speed of light, m the mass of binary, r and v the orbital radius and velocity; the characteristic radius for a compact 
object being R ps Gm/c 2 .) The post-Newtonian order is their ratio PNO(5f) = O e (8f)/O t (f). When the time derivative / is 
given, we can evaluate the order of the change Sf generated during a characteristic time-interval of the quasi-periodic motion: 
PNO(Sf) = O e (f)Oe(T)/O e (f). Here T is the radial period and O e (T) = re~ 1/2 /c. For a dimensionless quantity, typically an 
angle variable a, the simple estimate PNO(Sa) = O e (a)O e (T) holds. 

A numerical example shows that for a binary consisting of massive black holes with mass lOOM© and orbital radius of 10 4 
km, i.e. , one-hundredth of the Hulse- Taylor system, the post-Newtonian parameter is e w Gm/rc 2 ~ 10~ 2 . 

We also assume rapid rotation, thus the characteristic velocity of rotation V is commesurable with the velocity of light 
0(V/c) — 1. If this assumption is lifted, the spin-spin effects appear at higher post-Newtonian orders than is assumed here. 
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their directions, magnitudes and the angles, all enlisted in Table 1 will be needed to a precision at which precessional 
effects do not contribute. 

All spin terms in the radiation reaction up to the 3/2 post-Newtonian order were evaluated in ]l2| |, starting from 
the expressions of the radiated power and total angular momentum loss derived by Kidder |llf| and from the Burke- 
Thorne potential The employed averagin g m ethod, based on a suitably introduced radial parameter \ an( l on 
the residue theorem was described in detail in |L4| ■ 

The second-order spin contribution to the radiation reaction was considered recently in |l5| ] . There we have devel- 
oped the toolchest for computing secular radiative effects of spin-spin origin, essentially by proposing a description in 
terms of constants of the motion and the angular average L of the magnitude of the orbital angular momentum L(x)- 
The latter depends on the Keplerian true anomaly x through a spin-spin term: 



L( X ) = L- 



G[i 2 



■S1S2 sinKi sinK2{2Acos[x+ 2(00 — ip)] 



2c 2 L 3 
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The quantity A — (G 2 m 2 /i 2 + 2£ , L 2 //i) 1 / 2 is the magnitude of the Laplace- Runge-Lenz vector for a Keplerian motion 
characterized by the energy E and magnitude of orbital angular momentum L. The angles Ki — cos _1 (Si • L), (£=1,2) 
together with 7 = cos _1 (Si • S2). characterize the relative orientation of the angular momentum vectors Si, S 2 and 
L. 

This description enabled us to introduce the generalized true anomaly parametrization r(x) for the radial motion: 
(Eq. (3.7) of Ref. E5[). By employing x as radial variable, the average of a generic function / over a radial period T 



T 



r dx 



(1.3) 



is computed by Taylor expanding f 2 taken from the radial equation, by substituting 



dr 
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(r mlI being the turning points) and finally, by switching to the complex parameter z = exp(z%). Then the residue 



theorem is employed in the computation of the integral (1.3), which in the majority of cases is given by the residues 
in the origin E3 . 

As application, we have computed the secular losses of the energy E and of L. They depend on E, L, on the 
angles "0i subtended by the node line with the projections of the spins in the plane of the orbit and on the angles Ki 
and 7. The enlisted five angles are constrained by two algebraic relations |12|: 



= Si sin Ki cos -01 + S2 sin K2 cos ip2 
cos 7 = cos Ki cos K2 + sin k\ sin K2 cos , 



(1.5) 
(1.6) 



where A?/> = ip2 — 01- We also introduce here the notation ib = (0i + "02)/2. In order to have a closed system of 
differential equations, beside dE/dt and dL/dt computed injO the radiative change of the angles Ki and 7 is required. 
This was given to 3/2 post-Newtonian order accuracy in |12[|. It is the purpose of the present paper to compute the 
radiative evolution of the angles Ki and 7 up to the second post-Newtonian order by the inclusion of the spin-spin 
terms. 



TABLE I. The effect of precessions on angular momenta. 



quantity (/) 



leading order change O(Sf) origin of the leading order change constant at order 



ATT 



£ 

£ 3/2 
c 3/2 



Si , Si = Si/5, , 7 = cos- 1 (S 1 -S 2 ) 
L , L = L/L 

Ki = cos _1 (Si ■ L) 
L 

S, 



e 

£ 3/2 
,3/2 



spin-orbit 
spin-orbit 
spin-spin 
spin-spin 
higher order effects 
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We complete the task by passing through the following steps. A well-known result states |T(| that no leading-order 
gravitational radiation leaves an axisymmetric body, as its quadrupole moment is a constant. It is also known that 
the radiative loss dSi/dt in the spins of the axisymmetric neutron stars or black holes forming a binary is two post- 
Newtonian orders higher that the leading radiative loss dL/dt of the orbital angular momentum jnj]. In Sec. 2 we 
derive a related generic result, which holds for binaries with axisymmetric constituents. We show that axisymmetry 
implies that the leading order secular contribution to the radiative change of the spins (which would appear at the 
second post-Newtonian order after the leading radiation terms) averages out during one cycle of the quasi-periodic 
motion. This result simplifies the study of the secular evolution of the angles. Indeed we will drop all terms of the 
form dSi/dt multiplied with any quantity which is a constant to the required accuracy. 

The rest of the paper is organized as follows. In Sec. 3 we discuss the radiative change in the angles Ki and 
7 by considering the Newtonian (N), first post-Newtonian (PN), second post- Newtonian (2PN), spin-orbit (SO), 
spin-spin (SS) and leading order tail contributions. We show that the N, PN, 2PN and tail terms drop out as a 
consequence of the property of the respective parts in dJ / dt of being aligned with the Newtonian part of the orbital 
angular momentum. The SO terms were given previously in jL2j. We complete the list by computing in Sec. 4 the 
averaged SS contributions. They contain both two-body and self-interaction terms. 

As a result we obtain the system of equations which govern the radiative secular evolution of the relative orientations 
of the angular momentum vectors Si, S 2 and L up to the second post-Newtonian order. 



II. RADIATIVE LOSS IN THE SPINS 

In jl^] we have derived the radiative loss in the spin Si of the i th axisymmetric body, following |l7| . The integral 
over the volume of the body of the moment of the reaction force (the sign swapped gradient of the Burke-Thornc 
potential) gave: 

,/,Sii " 2G ^-Ae^I^SSM^- (2-1) 



dt ~ 5c^ J 

(There is no summation over i.) Each body is characterized by its principal moments of inertia 8^ and 0^ and by the 

angular velocity fl; = Si/Q'i, while iff 1 " 7 is the 5 th time derivative of the s ystem 's Newtonian symmetric trace- free 
mass quadrupole moment. As the loss in the spin vector described by Eq. ( |2.l[ ) is o f second post-Newtonian order 



(above the leading radiation term in the orbital angular momentum loss), Eq. (2.1) is still valid to this order for 
approximately axisymmetric bodies, with the deviation from axisymmetry being of any post-Newtonian order. We 
emphasize that Eq. ( pTl| ) implies dSi/dt = SidSi/dt, therefore to this order the radiation reaction will change the 
orientation but not the magnitude of the spin vectors. We also stress that dS-Jdt generates 3/2 post-Newtonian order 
radiative changes S ra d Si above the leading order radiation losses (in L for example). 

In coordinates (x,y,z) = r(cos - 0,sin'!/>,O) the spins are expressed as Si = Sj(sin k% cos ipi, sin Ki smipi, cos Ki). By 
employing the Keplerian equation of motion ajy = —Gmr/r 3 and the radial equation r 2 = 2E/ fi + 2Gm/r — L 2 / /j, 2 r 2 



in Eq. (2.1), we obtain the instantaneous spin-loss equation: 



1 d(S0„ 2G 2 msin K4 ( Q± y 



Si dt 5c \&i J "* 

Aii = cosKi[ai sin(2-0 — ipi) + °2 cos(2ip — ipi) + 03 sini/'t] 
Ai2 = cosKi[a2 sin(2-0 — ipi) — ai cos(2-0 — ipi) — 03 cos^i] 

Aj3 = - sinKj[ai sin2(V> - ipi) + ai cos2(i/> - ipi)] , (2.2) 



with the coefficients 



ai = firr(12fiEr 2 + 20Gm^i 2 r + 45L 2 ) 
a 2 = 4L(18/i£r 2 + 20Gm^ 2 r - 15L 2 ) 

a 3 = nrr(l2{iEr 2 + 20Gmfi 2 r - 15L 2 ) . (2.3) 

In order to obtain the secular radiative changes of the spins we insert the Newtonian expressions 

A L 2 

f = T sinx, r = —r^ r- 7 ip = X + ipo- (2.4) 

L /x(Gm/i + Acosx) 
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in Eq. Q. 

Then we average by the method described in the Introduction (see also fl5j] ) employing the Newtonian expression 
(l/r){dr / d\) = dt/dx = ^r 2 /L. We obtain the generic result that the radiative change in the spins averages to zero 
in the second post- Newtonian approximation: 

(§>->■ 

Therefore the averaging-out property of some particular projections of the spin losses found in jl2j and |l^] also 
emerge. 

As the first correction to both the acceleration and the Burke-Thorne potential is one post-Newtonian order higher 
than the respective leading terms Jl9| , the averaging-out property of the radiative change in the spins derived above 
holds at 5/2 post-Newtonian orders as well. 



III. RADIATIVE EVOLUTION OF THE ANGLES m AND 7 

We obtain the equations for the radiative evolution of ki = cos _1 (Sj • L), (i — 1,2) and 7 = cos _1 (Si • S 2 ) by 
differentiating the defining relations of these angles. 

The simplest is the evolution of the angle 7 spanned by the two spin vectors: 



-cos 7 =-(s a .S 2 ) = S 1 .— + S 2 — 



dS 1 



(3.1) 



The terms of the right hand side were already evaluated in the framework of the spin-orbit contributions in JL2| . 
Though they do not contain the spin explicitly, the angular velocities fii are present in their expression, Eq. (3.13) 
of [l2] ]. An order of magnitude estimate shows that the ratio of the usual spin-orbit terms (taken for example from 
the losses of Eqs. (3.10)-(3.12) of |12|) and such terms is of the order unity for rapidly rotating objects. It was 
computed in p"2[ , and one can check by simple inspection of Eq. (2.5) that 7 will receive no secular radiative change. 
The angles spanned by the spins with the orbital angular momentum evolve in a more complicated fashion: 



d _ d /- L \ _ 1 

dt C ° SKi -dt{ bi -W))~L(x) 



- d(J - S x - S 2 ) dL _ dSi 

Si-— r — ~ -r cos Ki + L~ 

dt dt dt 



(3.2) 



In the forthcoming expressions we replace L(x) L m all post-Newtonian terms. This is possible in the required 
second order accuracy as L(x) differs from the constant L in spin-spin terms (see Eq. (1.2)). 

By inserting the following expression for the loss of the magnitude of orbital angular momentum 



dL _^ <iL _^ dJ 

dt dt dt 



L 



dS 1 
~~dT 



dt 



(3.3) 



into Eq. (|3.2|), and from Si • dSi/dt = (no summation) we obtain: 



d 1 ~ dJ 1 

— COSKi = ——(Si - LcOSKj)-— + T 

dt MX) dt L 



~ rfSj ~ / dSi dS 2 \ , T rfS ; 



(3.4) 



where j ^ i. Remarkably all terms in the square bracket average out due to Eq^^). Therefore we will not consider 
them further. 

For the evaluation of the first term of Eq. ( |3.4[ ) we need the loss in the total angular momentum dJ/dt. It has the 
following structureR O] : 



2 To avoid confusion in notation, we denote the post-Newtonian, 2 d post-Newtonian, spin-orbit and spin-spin terms from a 
decomposition of a quantity expressed in terms of TjV, S; and r by lower-case indices pn, 2pn, so and ss, respectively. This 



notation applies to all terms of dJ/dt computed in |idj and JtJ . After inserting the expressions for v 2 , r 2 and 
(2.33) of [[L5j and (2.4), the terms of the new decomposition will be denoted by the respective capital letters. 



Eqs. (2.32), 
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dJ 

~dt 



N 



2pn 



(3.5) 



tail 



The leading order contribution was computed long time ago by Peters [|3| , the pn contribution by Junker and Schafer 
p2fl . the so contribution by Kidder and the 2pn contribution by Gopakumar and Iyer |t]] in terms of r, v, S ; and 
f. The method of computing the ss contribution was also indicated in |1 1| . The leading tail contribution appears at 
3/2 post-Newtonian order and its time average was given by Rieth and Schafer |2(]] in the form of a Fourier series. 
A property we want to stress is the following: 



~dt) 



= r , i, 2 -Ljv (3.6) 

N, pn, 2pn 

where the coefficients Tq i 2 can be read from Eq. (3.28.a,b) of and from (3.9.d) of §. The same property applies 
to the averaged tail contribution, given by Eq. (84) of |2Cj . To see this we note that all coefficients (denoted there by 
n (So)) of the spectral decomposition of the angular momentum loss (denoted there by (S)) have only z— components, 
thus: 



dJ 

dt 



tail 



L 



JY 



tail 



The next essential remark is that from 

L = L^r + hpN + L2PN 4 
we can express the orbital angular momentum as 

Ljv = (1 - 7i 



Lso = (1 + 7i + 72)Lat + L so 



72 + 7i)L - L 



so 



(3.7) 



(3.8) 



(3.9) 



Here the coefficients 7i ] 2 are given by Eqs. (2.9.b,d) of [[H] and they are of first and second post-Newtonian order, 
respectively. Therefore up to an accuracy of second post-Newtonian order we can write 

= (To+m + ?/2)L-r • Lc 

iV+ pn + 2pn 



dt 



■'SO 



(3.10) 



and 



tail 



tail 



where the coefficients of the expansion (3.10) are given by 

r)i = Ti - T 7i 



m = r 2 - Ti7i - r (72 - 7?) 



(3.11) 



(3.12) 



We w ill not need their explicit expressions. By inserting Eq. ( 3.10 ) in Eq. (3.5), then the resulting expression in Eq. 
(|3.4[), due to the remark 



(Si — LcosKi) ■ L =0 



we obtain for the first term of Eq. (3.4): 



(Si — L COS Hi) 



= ^(Si — LcosKi) 
L 



tail 



(3.13) 



(3.14) 



The tail term averages out due to Eq. ( 3.11 ) and property ( 3.13 ). 

As suggested by the notation, the second and third terms of Eq. (3.14) are those spin-orbit contributions to the 
radiative loss of the angles k^, which do not originate in the Burkc-Thorne potential. Indeed, modulo Burke-Thorne 
type terms they can be put into the concise form (Si • dL/dt)so which was computed previously (Eqs. (3.8)-(3.10) 
of |lj]). The averaged expression for the spin-orbit type loss in k,\ was also given by Eq. (4.4) of [Q, and a similar 
expression can be found for th e loss of K2 by interchanging the indices 1 <-> 2 and the ratios 77 = mijra\ «-> 17 . 

The last term of Eq. (3.14), modulo Burke-Thorne type terms, is the spin-spin part of the radiative evolution of 
the angles Kf. 



d 

— COS K, 

dt 



L cos/,-,- i • ( ^ 



(3.15) 



which will be computed in the next section. (We have denoted by ~ the equality modulo Burkc-Thorne type terms.) 
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IV. SECULAR RADIATIVE EVOLUTION EQUATIONS 



We start from the expression of (dJ/dt) ss given by Kidder JliJ in terms of the time derivatives of the mass quadrupole 
and velocity quadrupole moments (see also Eq. (4.18) of JTa] for the respective expression with c ^ 1 =/= G). The 
required SO part of the velocity quadrupole moment was computed first by Kidder and later verified by several 
authors. (Rieth and Schafer p0| presented a derivation based on a different spin supplementary condition while Owen, 
Tagoshi and Ohasha pjj have employed a ^-function type energy-momentum tensor.) 

After computing the expression (dJ / dt) ss in detail, we rewrite the last term of Eq. ( 3.14 ) as function of the radial 
variables r(\) and X- The procedure we follow was described in detail in fl5|| . As a result both self-interaction and 
two-body spin-spin terms emerge: 



-=- (Si — L cos Ki) 



1 ,A 



2G 3 m 2 /u sin m 



SS — self 



S1S2 



5c 7 r 6 
2G 2 S"iS* 2 sinK 



mi 



sm Ki cos Ki 



sin Kj cos Kj cos Aip 



(4.1) 



1 Ki r 

Tj — |ui[sm re. 



5c 7 fj, 2 L 2 r 
-U2 [sin Ki cos Kj cos 2(\ + ipo — ipi) 
- sin x{ u 3 sm K j cos K i sin-iipj ~ ip%) 
-U4[sin Ki cos Kj sin 2{% + ipo — ipi) 



cos Kj + sin Kj cos cos Aip] 

sin Kj cos Ki cos 2 (% 



sin Kj cos Ki sin 2(% 



- ip - ip)} 
iPo~4>)]}} ■ (4-2) 



Here j ^ i and the coefficients are 



ui = -L 2 (\2^Er 2 + AGmn 2 r - 15L 2 ) 

u 3 = 3fiAr{~2^Er 2 + 5L 2 ) 
Ui = 3(iAr(-2i_iEr 2 + 3L 2 ) 



(4.3) 



The averaging procedure based on the parametrization r{\) and on the residue theorem yields the self-interaction 
and two-body spin-spin terms in the secular radiative loss of the angles Ki and 7. We enlist them together with the 
spin-orbit contributions: 



=+) =0 



d~) 
~dt 

dre, 

~~dt 



dKi 
~dt 



d^ 
~dt 



dKi 
~~dt 



SO \ / SS-self 

given by Eq. (4.4) of [12] 



dKi 
~~dt 



ss 



so 



I dKi 
\~dt 



SS — self 



\~dt 



G 2 ra V i{-2 [ iE fl 2 1 
20c 7 Z 9 

G 2 (-2 M £;) 3 / 2 5 1 5 2 



sm Ki cos re» 



sin cos Kj cos A-0 



Si5 2 



I V2 (sin Ki cos Kj + sin Kj cos Kj cos A?/ 1 ) 



20c 7 L 9 sin re, 
+ V3 [sin Ki cos Kj cos 2 (-00 — ip 
where j ^ i and the coefficients V1-3 are: 

60G 2 m V#Z 2 + 35G 4 mV 



sin Kj cos Ki cos2(t/>o — VO]} 



(4.4) 
(4.5) 
(4.6) 
(4.7) 

(4.8) 



Vi = 12E 2 L 4 -+ 
V 2 = 564£ 2 Z 4 + 1620G 2 m 2 (i 3 EL 2 + 805G 4 mV 
V 3 = 60{8E 2 L 4 + l8G 2 m 2 n 3 EL 2 + 7G 4 m 4 u 6 ) . 



(4.9) 

Note that Vi = -Di, because both arose from the average of r -6 . (The coefficient D±, given by Eq. (4.32) of p5[ , 
governs the self-interaction term in the radiative loss of L). In order the 3/2 post-Newtonian order- accurate spin-orbit 
contribution to hold at 2PN, the replacements L —> L and Aq — > ^4 should be carried on in Eq. (4.4) of H . 

We stress that Eqs. (4.4)-( |Q| ) contain radiative terms in the angular evolution up to the second post-Newtonian 
order above the leading radiative effects. 
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V. CONCLUDING REMARKS 



By proving the 5/2 post-Newtonian accurate result that axisymmetric objects do not radiate away any fraction of 
their initial spins, the computation of the secular angular evolutions induced by radiation reaction became simpler. 
Then from the analysis of dJ/dt we could conclude that there are no N, PN, 2PN and tail contributions to the 
secular radiative angular evolutions. The SO part of the equations was given in |Q . We have derived the SS terms 
here. 

No spin-spin terms appear in the radiative evolution of the angle 7. As the spin-orbit terms average out p2| , 
we found the remarkable result that the angle spanned by the spins receives no radiative secular change up to the 
second post-Newtonian order. (However the angle 7, together with all other angles is subjected to precessional - 
both instantaneous and secular - evolution. Although the second post-Newtonian order precessional (nonradiative) 
evolution of the angles is not developed in the paper, inspection of Table 1 shows that the precessions do not contribute 
to the 2PN radiative evolution of the angles /t; and 7, which contain only 3/2 PN and 2PN parts.) 

The angl es K j evo lve under the influence of both spin-orbit and s pin- spin terms. The radiative angular evolution 
equations (|4.4| )-(4.g ) together with the algebraic constraints (L5)-([Tli|) and the expressions for dE/dt and dL/dt 
derived in |15| form a closed system of first order differential equations. 

The angle tpo appearing in this system is an integration constant. It is interpreted as the angle subtended by the 
node line with r at \ — (with the periastron line). Each precession modifies the value of ipo by a small amount. 
According to the arguments of Ryan in pq] , where this angle first appeared, the terms containing periodic functions 
of ipo average to zero whenever the precession time scale is short compared to the radiation reaction time scale. 

As it happened with the loss of E and L Jl5| , the spin-spin terms of the radiative Ki-evolution could be decomposed 
into two-body and self-interaction terms. In the one-spin limit S2 = the terms proportional to Sf from Eq.(1.7) 
represent the second post-Newtonian correction to the radiative evolution of the angle n\ derived earlier in the 
Lense-Thirring approximation fl24| . 
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